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Extending to continuous potentials a cleaving wall molecular-dynamics simulation method recently 
developed for the hard-sphere system [Phys. Rev. Lett 85, 4751 (2000)], we calculate the crystal-melt 
interfacial free energies, 7, for a Lennard-Jones system as functions of both crystal orientation and 
temperature. At the triple point, T* = 0.617, the results are consistent with an earlier cleaving 
potential calculation by Broughton and Gilmer [J. Chem. Phys. 84, 5759 (1986)], however, the 
greater precision of the current calculation allows us to accurately determine the anisotropy of 7. 
From our data we find that, at all temperatures studied, 7111 < 7110 < 7100- Comparison is made 
to the results from our previous hard-sphere calculation and to recent results for Ni by Asta, Hoyt 
and Karma [Phys. Rev. B, 66 lOOlOl(R) (2002)]. 



INTRODUCTION 

The magnitude and orientational dependence 
(anisotropy) of the solid-liquid interfacial free energy, 
7si, is a primary controlling parameter in the kinetics 
and morpholgy of crystal growth from the melt[? ], 
especially in the case of dendritic growth [? ]. As a 
consequence, the ability to accurately measure or pre- 
dict this quantity for specific materials is of significant 
technological and scientific importance. For most ma- 
terials, the only experimental data for 7^; is extracted 
indirectly from nucleation data (assuming some level 
of classical nucleation theory) [? ? ? ]. Such indirect 
measurements tend to underestimate the actual inter- 
facial free energy by 10-20% and represent orientatial 
averages, so all information as to interfacial anisotropy 
is lost. Direct experimental measurements, usually in- 
volving contact angle studies, are quite difficult and 
relatively few in number[? ] and, with the exception 
of a small number of studies on transparent organic ma- 
terials [? ? ], are not of sufficient precision to resolve 
anisotropy. This paucity of reliable direct experimen- 
tal measurements on technologically useful materials 
(such as metals) has motivated the development of a 
variety of novel computational methods to determine 
7si via molecular simulation[? ? ? ? ]. 

The interfacial free energy of a crystal-melt interface 
is defined [? ] as the reversible work required to form 
a unit area of interface. In a simulation this can be 
accomplished by constructing a continuous thermody- 
namic path between an initial system consisting of sep- 
arated bulk crystal and liquid to a final state contain- 
ing an interface. The value of jsi is then determined by 
thermodynamic integration[? ] along that path. This 



is a tedious process and care must to be taken to ensure 
that the process is reversible, i.e., intergration along 
the path in both the forward and backward directions 
yields the same result (no hysteresis). The first such 
calculation on a crystal-melt interface was performed 
by Broughton and Gilmer [? ] on a system of parti- 
cles interacting with a LJ potential (truncated so that 
both the potential and the force vanish at 2.5(7, where 
a is the usual LJ diameter). To perform the thermody- 
namic integration they employed external cleaving po- 
tentials that were slowly turned on to separate the sam- 
ples. The precise (rather complex) functional forms of 
the cleaving potentials were chosen, more or less by 
trial and error, to minimize hysteresis. The values of 
jsi were determined to be 0.35(2), 0.34(2) and 0.36(2) 
(in units of e/cr^) for the [111], [100], and [110] crystal 
orientations, respectively. The numbers in parentheses 
give the uncertainties in the last digit shown. The main 
source of this error is the small amount of hysteresis in 
the cleaving of the liquid phase. Unfortunatly, the pre- 
cision of these results was not sufficient to resolve the 
anisotropy of the interfacial free energy for this system. 

Recently, we have developed a modification of the 
Broughton and Gilmer approach in which planar cleav- 
ing walls, as opposed to cleaving potentials, are used to 
separate the phases. These walls are constructed out 
of the same type of particles as present in the system, 
with a 2-d geometry consistent with the symmetry of 
the interfacial orientation under study. This method 
was originally applied to the crystal-melt interface of 
a system of hard spheres [? ] and was shown to have 
sufficient precision to resolve the anisotropy with [111] 
being the lowest free energy face, followed by [100] and 
[110]. The cleaving wall method is complementary to 



a method due to Hoyt, Asta and Karma[? ] in which 
the interfacial free energy is extracted from the interfa- 
cial stiffness, which is determined from the spectrum of 
fluctutations in interfacial position. The cleaving-wall 
method has an advantage in that it requires simula- 
tions with an order of magnitude fewer particles than 
those required for the fluctuation method (^^[lO*] vs. 
©[10^]). The precision in the raw values of the inter- 
facial free energies is greater for the the cleaving wall 
method than in the fluctuation approach (even con- 
sidering the far smaller system samples); however, the 
fluctuation method yield somewhat more precise values 
of the anisotropy parameters since these are obtained 
from the anisotropy of the interfacial stiffness, which is 
far more anisotropic than the interfacial free energy. 

In this work, we extend our cleaving wall approach to 
systems of particles interacting with continuous poten- 
tials, specifically applying it to the system of truncated 
Lennard-Jones (LJ) particles considered by Broughton 
and Gilmer. Our results at the triple point are con- 
sistent with Broughton and Gilmer's calculation, but 
are of greater precision allowing us to resolve the 
anisotropy, which is found to differ slightly from that 
for the hard-sphere system[? ]. In addition, we de- 
termine the temperature dependence of jsi along the 
coexistence curve. The magnitude of jsi is shown to 
scale roughly linearly with the melting temperature (as 
predicted by the hard-sphere model[? ]). 

THE CLEAVING WALL METHOD 

The direct determination of the excess free energy 
of the crystal-fluid interface of a model system within 
a computer simulation can be achieved by thermody- 
namic integration along a reversible path beginning 
with separate crystal and fluid bulk systems prepared 
at the coexistence conditions and ending with a sys- 
tem containing a crystal-liquid interface at equilib- 
rium with the surrounding bulk phases. Construction 
of such a path requires the development of a proce- 
dure to reversibly cleave a simulation box into two 
non-interacting systems. Following the prescription of 
Broughton and Gilmer [? ], we identify the following 
steps in the process of creating the crystal-liquid inter- 
face (see Fig. P): 

Step 1: Split the crystal bulk system with a suitably 
chosen "cleaving" potential while maintaing the 
periodic boundary conditions 

Step 2: Split the liquid system in a similar way 

Step 3: Juxtapose the cleaved crystal and liquid sys- 
tems by rearranging the boundary conditions 
while maintaining the cleaving potentials 
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FIG. 1: Illustration of the four steps reversible process 
of creating the crystal- liquid interface from separate bulk 
systems. Dashed lines show the location of the cleaving 
planes and thick solid lines represent the cleaving potental. 
Letters A and B label the crystal and liquid on either side 
of the cleaving planes for better illustration of the bound- 
ary conditions rearrangement in step 3. Periodic boundary 
conditions are assumed on all boundaries of the simulation 
boxes. 



Step 4: Remove the cleaving potentials from the com- 
bined system. 

The interfacial free energy is calculated as the total 
work required to perform the above steps divided by 
the area of the constructed interface. In addition to 
the coexistence conditions, the result will also depend 
on the orientation of the crystal with respect to the 
interfacial plane. We will refer to the plane along which 
the crystal and liquid systems are split as the cleaving 
plane. The location of the cleaving plane in the crystal 
system is chosen in the center of the simulation box 
between two crystal layers, while in the liquid system 
the precise location is arbitrary. 

The reversibility and precision of the thermodynam- 
ics integration process are very sensitive to the choice 
of cleaving potentials. The requirements for the cleav- 
ing potential are two-fold: First, the cleaving potential 
should perturb the system as little as possible. As a 
consequence, it is desirable that, in step 2, the potential 
introduces structure into the cleaved liquid that is com- 
patible with the structure of the crystal layers. Sec- 
ond, the cleaving potential must be strong enough to 
prevent the particles from crossing the cleaving plane. 
Otherwise, the rearrangement of the boundary condi- 
tions in step 3 cannot be performed. 

Broughton and Gilmer [? ] carefully designed a set 
of cleaving potentials for the truncated Lennard-Jones 
(LJ) system. However, their approach for constructing 
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the cleaving potential was optimized specifically for the 

LJ system and is not easily adaptable to a general case. 
Here we outline an approach that is easily adaptable 
to systems with different interparticle interaction po- 
tentials and different crystal structures. 

As was demonstrated in our recent calculation of the 
interfacial free energy for the hard-sphere system[? ], 
the appropriate structure in the interfacial region of the 
fluid can be easily introduced by the potential of a pair 
of "cleaving walls" , made of properly oriented crystal 
layers, each interacting only with the particles on the 
opposite side of the cleaving plane. The wall particles 
are held fixed at the crystal lattice sites. When the 
two walls arc far away from the cleaving plane, they do 
not interact with the system. Moving the walls closer 
to the cleaving plane introduces a cleaving potential 
whose structure is similar to that outside the crystal 
layers of the studied system. When a liquid system 
interacts with such a cleaving wall, it is expected to 
form an interfacial layer structure similar to that at the 
crystal-liquid interface. To achieve this, the interaction 
potential of the wall particles must be similar to that 
of the system particles. 

Since interactions in the crystal are usually domi- 
nated by the short-range repulsive part of the potential, 
it is sufficient to choose the interaction potential of the 
wall particles as a monotonically decreasing function 
(l){r) with a relatively small cut-off radius r^,, which 
models the repulsive core of the interparticle potential 
in the system under investigation. Interaction of the 
system with each cleaving wall is then given by 

$i(r;z)=^</.(|r-R«+nz|), 
j 

$2(r;^) = 5^<^(|r-Rf -n^l), (1) 
j 

where R-j^'^^ are the positions of the particles forming 

the walls, n is a unit vector normal to the cleaving 
plane, and z measures the distance of the walls to the 
cleaving plane. 

Next, we need to ensure that the system on the one 
side of the cleaving plane interacts only with the wall 
on the other side of the plane. To achieve this, we 
use the monotonic character of the potential (f>(r) and 
define the cleaving potential as the minimum of the 
two wall potentials, namely 

$(r;z) = min($i,$2), (2) 

which decays to zero away from the cleaving plane, 
as long as the potential (j){r) is chosen to be a pos- 
itive monotonically decreasing function. To remove 
discontinuity of the gradient of $ at the points where 



^1 = ^2, we slightly modify the minimum function as 
follows: 

{X, x<y-S{x + y) 

y, x>y + S{x + y) (3) 
p{x,y), otherwise 

with 

and parameter 6 characterising the relative width of the 

interpolation region. We set S = 0.25 for the present 
study. The cleaving potential is defined as 

*(r;z) =rn($i,$2)- (5) 

Even though the crystal system does not require ad- 
ditional ordering and, as proposed by Broughton and 
Gilmer, can be cleaved with a short-range repulsive po- 
tential centered at the cleaving plane, we have found 
that using the same cleaving potential for both crystal 
and liquid systems results minimizes the error during 
the thermodynamic integration in step 3. 

To calculate the reversible work in steps 1, 2, and 
4, we can use the wall position z as the integration 
coordinate. The reversible work is thus determined by 
evaluating the integral 

where the angle brackets denote averaging over a sim- 
ulation run at a fixed cleaving wall position. In steps 
1 and 2, the initial position of the cleaving walls, Zi, 
is just outside the range of the interaction potential 
determined by the cut-off radius r^. The final wall 
position Zf is determined by the requirement that the 
cleaving potential is sufficiently strong to prevent the 
particles from crossing the cleaving plane. In step 4, 
the initial and final positions of the walls are reversed. 
Because of the repulsive character of the cleaving po- 
tential, the work in steps 1 and 2 is expected to be 
positive, while in step 4 it is negative. 

In step 3, the boundary conditions are gradually re- 
arranged using a coupling parameter A. The total in- 
teraction energy in step 3 is given by 

U{X) = (1-A) Yl ^(^ij) 

i<j, AB 

+A + (7) 

i<j, AA i 

where the letters (AB) refer to the boundary con- 
ditions with crystal- crystal and liquid-liquid interac- 
tions across the cleaving planes, and {AA) refer to 
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the boundary conditions with crystal-Uquid interac- 
tions (see the diagram). The last term represents the 
cleaving potential at the final wall positions in steps 
1 and 2. The work done during step 3 is determined 
from the integral 



(8) 



where 



dX 



i<j, AA i<j, AB 



RESULTS FOR TRUNCATED 
LENNARD-JONES POTENTIAL 

In order to have direct comparison with the results 
of Broughton and Gilmer[? ], we have used the same 
modification of the Lennard-Jones potential, namely [? 



u{r) = 



4e 
0, 



r 

, 12 



Ci r< 2.3a 

2.3(7 <r < 2.5a 
2.5a < r. 



where Ci = 0.016132e, C2 = 3136.6e, C3 = -68.069e, 
C4 = -0.083312e, and C5 = 0.74689e. (Note: the sign 
of C4 was incorrectly reported as positive in Broughton 
and Gilmer's original publication). This potential has 
a continuous first derivative and a small discontinuity 
in its second derivative at r = 2.3a. 

To achieve optimal performance of the simulation at 
all stages of the cleaving process, we use a cell-assisted 
force evaluation method [? ]. The simulation region 
is subdivided into cells with edge sizes just exceed- 
ing the interatomic interaction range. When atoms 
are assigned to cells according to their current posi- 
tion, only interactions between atoms in the same cell 
or in immediately adjacent cells needs to be consid- 
ered. Such a method is particularly useful in step 3, 
where interaction forces must be computed for two dif- 
ferent boundary conditions [see Eq. (TJ]. With the cell- 
assisted method the additional computational effort is 
limited to the cell layers immediately adjacent to the 
cleaving planes. (Note, the cell subdivision is carried 
out in such a way that the cell boundary coincides with 
the cleaving plane.) 

The insertion of the exteral cleaving potentials 
("walls") breaks the translational symmetry of the 
system Hamiltonian, linear momentum is no longer 
strictly conserved. This leads to a problem for long 



runs in that the crystal, taken as a whole, can drift rel- 
ative to the cleaving plane. To prevent this, we immo- 
bilize the innermost 2 crystal layers by assigning them 
an infinite mass. In the data-collection runs these fixed 
planes were 20 crystal layers away from the nearest 
crystal-melt interface, and should not have a signifi- 
cant effect on the interfacial properties. To check this, 
we have performed smaller simulations with identical 
cross-sectional area, but fewer crystal layers (where the 
interface is about 12-14 crystal layers distant from the 
fixed planes) and obtain results identical, within simu- 
lation error, to the larger systems. 

For the present study, we calculate 7 at both the 
triple point temperature of T* = ksT/e = 0.617 (as 
determined by Broughton and Gilmer [? ]), and at 
higher temperatures {ksT/e = 1.0 and 1.5). The 
crystal-liquid coexistence parameters at these temper- 
atures are summarized in Table U (Note the slightly 
negative pressure at ksT/e = 0.617. We have found 
that kBT/e = 0.618 is a better estimate of the triple 
point temperature for the modified LJ potential used 
here. However, to have a direct comparison with the 
results of [? ], we have calculated interfacial free en- 
ergy at the lower temperature.) The parameters are 
obtained by monitoring long simulation runs of the 
crystal-liquid interfacial systems. If the initial condi- 
tions are shifted from those at coexistence, such a sys- 
tem equlibrates itself through melting/freezing at the 
interfaces. During the runs, the pressure tensor profiles 
[? ] are monitored and the simulation box rescaled, if 
necessary, to remove any stress in the bulk crystal. 



TABLE I: Coexistence conditions for the modified 
Lennard-Jones potential 



r, fcse-i 


Pc, cr-^ 


Pi, 


P, eo-^ 


0.617 


0.945 


0.828 


-0.02 


1.0 


1.005 


0.923 


4.95 


1.5 


1.074 


1.003 


12.9 



At each of the three temperatures, we calculate the 
interfacial free energy for three crystal orientations: 
[100], [110], and [111]. The cleaving wall is constructed 
from a single crystal layer for the [100] and [111] orien- 
tations and from two layers for the [110] orientation 
(this was necessary to prevent atoms from crossing 
through this non-close-packed face). As discussed in 
the previous section, the cleaving potential 4'{r) is con- 
structed from the repulsive core of the LJ potential: 
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FIG. 2: Integrand (Eq. 6) for thermodynamic integration 
in steps 1,2 and 4 for the [111] interface at T* = 0.617. 
The error bars are two standard deviations plus estimated 
hysteresis error, where appropriate. 

The results are summarized in Tabled For illustration 
of the method we show in Fig [21 the thermodynamic 
integration integrand (Eq. 6) for steps 1,2 and 4 for the 
[111] interface at T* = 0.617. For that same interface 
and temperature, the integrand for step 3 (Eq. 7) is 
shown in Fig. |31 with an inset highlighting the region 
of maximum hysteresis error. 

TABLE 11: Interfacial free energy (in units of ecr^'^) for 
the truncated Lennard- Jones potential for selected temper- 
atures and crystal orientations. Numbers in parentheses in- 
dicate two standard deviation error bars (plus an estimate 
of the hysteresis error on the last digit (s) shown. 





T* = 0.617 


1.0 


1.5 


7100 


0.371(3) 


0.562(6) 


0.84(2) 


7110 


0.360(3) 


0.543(6) 


0.82(2) 


7111 


0.347(3) 


0.508(8) 


0.75(3) 



Even though the relative statistical accuracy that 
can be achieved for a given duration of the simulation 
run is about the same for all temperatures and orien- 
tations (approximately 0.5%), we see that the relative 
error range in Table ITU increases with the temperature, 
especially for the [111] orientation. This increase is 
due to the observation that reversibility (as measured 
by the lack of hysteresis) of the thermodynamic inte- 
gration process for the LJ system is more difficult to 
achieve for higher temperatures. 

There are two sources of potential hysteresis in the 
four-step process of creating an interface. The first one 
is associated with the liquid ordering transition that 
occurs either at the end of step 2, or at the beginning 
of step 3. For lower temperatures, the hysteresis can be 
essentially eliminated by increasing the duration of the 
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FIG. 3: Integrand for thermodynamic integration in step 
3 for the [111] interface at T* — 0.617. The inset shows a 
magnification of the region primarily affected by hysteresis, 
reflected by the larger than average error bars. The error 
bars are two standard deviations plus estimated hysteresis 
error, where appropriate. 

equilibration runs near the point of transition. How- 
ever, the hysteresis for higher temperatures is more 
persistent, especially in the case of the [111] orienta- 
tion. This may be due to a particular choice of the 
cleaving potential. More research is necessary to eluci- 
date the origin of the hysteresis and ways to eliminate 
it. 

The second source of hysteresis is the fluctuation in 
the interface position at the end of step 4. When the 
cleaving potential is removed from the interfacial sys- 
tem in step 4, the position of the interface is no longer 
tied to the cleaving plane. Because the system contains 
two interfaces, they can change their position without 
disturbing system equilibrium by the process of melting 
at one interface and simultaneous freezing at the other 
interface. The difficulty in verifying the reversibility of 
step 4 is that the mobility of the interfaces causes the 
reverse process to follow a slightly different thermody- 
namic integration path. To deal with this problem, we 
try the reverse step 4 process on the interfacial systems 
after several equlibration runs of various duration and 
select the one with the path closest to the forward pro- 
cess. The difference in the calculated work during the 
forward and reverse processes is accounted for in the 
estimate of the error range given in Table ^ 



ANALYSIS AND SUMMARY 

The error bars on the calculations described above 
are small enough to resolve the anisotropy in the inter- 
facial free energy for the Lennard- Jones system. While 
we have determined 7s; only for the [100], [110] and 
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[Ill] directions, it is possibile to extract from these 
data some information as to the full angular depen- 
dence of the free energy. Defining the orientation unit 
vector n as the unit vector perpendicular to the in- 
terfacial plance, one defines an orientation dependent 
interfacial free energy 7(n), which can be parameter- 
ized by an expansion in terms of cubic harmonics. One 
such expansion, due to Fehlner and Vosko, has been 
recently applied to the interfacial free energy of Ni/Cu 
alloys by Asta, et al. [? ] In terms of the Cartesian com- 
ponents of n = {ni, 712, 713} , this expansion (truncated 
at sixth-order) is 



7(n) 



70 




(11) 



-£2 



E 



4 , r-r; 2 2 2 



where 70 is the orientationlly averaged interfacial free 
energy and ei and €2 are expansion coefficients. This 
parameterization has an advantage over the so-called 
"Kubic Harmonic" expansion [? ] that has been used 
recently to parameterize free energy anisotropy in met- 
als [? ? ] in that the expansion terms are orthogonal[? 
] . For the orientations studied here this expansion be- 
comes 



7100 
7110 
7111 



2 4 . 
5^1 + 7^2] 

1 13 

10'^ " 14'^] 
n 4 64 . 
7o[l - T^ei + ^£2] 



70 [1 
70 [1 



15 



63 



(12) 



Using our values for 7 for the [100], [110] and [111] 
orientations one can solve for the three parameters 70, 
ei and £2: 



70 
ei 

€2 



IO7100 + I67110 + 97 



111 



35 

357100 - 87110 - 277111 
2270 

3 

(7100 - 47110 -I- 37111) 



(13) 



2270 



The anisotropy parameters for the Lennard-Jones 
system studied here, as well as those for the hard- 
sphere system calculated from the data from our ear- 
lier calculation [? ], are summarized in Table UTTl Also 

included in Table IIIII are the values of and 
— 70 

'^""70^"^ , which also serve to quantify the anisotropy. 

The error bars in 70 are smaller than those for the raw 

interfacial free energies since 70 represents a weighted 

average of similar numbers - a process that decreases 



statistical error. The anisotropy parameters repre- 
sent differential quantities involving diff'ferences be- 
tween similar numbers (which magnifies relative error) 
so the relative error bars for those quantities are larger 
than in the raw data. This is in contrast to the fluc- 
tuation method[? ] where the relative error in the in- 
terfacial free energy is larger than in the current study 
(despite the much larger systems), but the anisotropy 
parameters are more precisely determined due to the 
fact that in that method the anisotropy is determined 
from the anisotropy in the interfacial stiffness, which is 
far greater than that of the interfacial free energy itself. 
From the anisotropy parameters given in Tablc Hlll the 



TABLE III: Interfacial free energy anisotropy parameters 
for the truncated Lennard-Jones at three different temper- 
atures. For reference, the same quantities for the hard- 
sphere system [Reference 9] are also included. Numbers in 
parentheses reflect 2-sigma error bars in the last digit(s). 
The units for 70 are e/cr^ for the LJ data and kT/a^ for 
hard spheres. 





LJ (this work) 
T* = 0.617 1.0 1.5 


HS ^ 


70 

£2 

7100 —7110 
70 

7100 —7111 
70 


0.360(2) 0.539(4) 0.808(13) 
0.093(17) 0.13(3) 0.15(6) 
-0.011(4) -0.022(9) -0.03(2) 
0.03(1) 0.035(15) 0.025(3) 
0.07(1) 0.10(2) 0.11(3) 


0.617 
0.07(3) 
-0.044(12) 
-0.032(22) 
0.065(22) 



trend is observed that the anisotropy parameters in- 
crease in magnitude with increasing temperature, with 
ei becoming more positive and £2 becoming more neg- 
ative. In comparison with the hard-sphere values, the 
value of £2 approaches the hard-sphere value as T gets 
larger, but the value of £1, which is at T* = 0.617 quite 
close to the hard-sphere value, diverges away from the 
hard-sphere value as T increases. Of course, in the 
limit of high temperatures, the Lennard-Jones system 
approaches an inverse twelth-power repulsive poten- 
tial, not the hard-sphere potential, so one would ex- 
pect the high temperature behavior of the anisotropy 
to approach that of the former potential. This implies 
that any perturbation theory for the interfacial free en- 
ergy with a hard-sphere reference potential will not be 
adequate to predict anisotropy, and that a study of the 
anisotropy for a variety of different possible repulsive 
potentials would be very useful. 

In recent work[? ], we have pointed out that the 
interfacal free energy for simple systems with face- 
centered cubic (fee) crystal structures can be quanti- 
tavely described by a hard-sphere model. It is useful 
to check this hypothesis here, since the LJ system is a 
standard molecular model that freezes to an fee crystal. 
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and since the interfacial free energy was determined in 
this study directly and not indirectly from nucleation 
data. The hard-sphere model predicts 

joiHS) = 0.617fcr/cr2 . (14) 

where cr^ is the hard-sphere diameter. For our 
Lennard-Jones system we can define an effective (tem- 
perature dependent) hard-sphere diameter using the 
Barker-Henderson criterion [? ] from liquid-state per- 
turbation theory 

/•OO 

'^eff = / {1 - exp[-Ur{r)/kT]}dr , 
Jo 

where Ur{r) is the repulsive part of the potential, which 
we define in the Week-Chandler- Anderson sense [? ] as 
the full potential truncated (and shifted) to zero be- 
yond the minimum of the attractive well. This pro- 
ceedure yields values of tTe// of 1.032cr, l.OlGcr and 
l.OOOcr, for T* = 0.617, 1.0 and 1.5, respectively In- 
serting these values into Eq. El yields predicted val- 
ues of 70 (in units of e/cr^) of 0.36, 0.60 and 0.93, for 
T* = 0.617, 1.0 and 1.5, respectively. The agreement 
with the values listed for LJ in Table IIIII is excellent 
at the lower temperature, but overestimates the ac- 
tual value by several percent at the higher tempera- 



tures. This agreement gives more evidence to support 
the general hypothesis [? ] that the interfacial free en- 
ergy of close packed systems is largely determined by 
packing considerations, not energy. 

It is interesting to note that the anisotropy param- 
eters ei and €2 for the truncated LJ potential at the 
triple point (T* = 0.617 and essentially zero pressure) 
are identical to those calculated[? ] for Ni at 1.00 atm 
(also essentially zero pressure). The fact that the pa- 
rameters are exactly the same for these two (essentially) 
zero pressure systems is, given the error bars, proba- 
bly coincidental; however, the data does show that the 
anisotropy for Ni is better modelled by a LJ poten- 
tial than by a hard-sphere potential. More study is 
required to determine the exact role that details of the 
potential play in determining interfacial anisotropy. 
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